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Abstract 

In this study, using the concept of relative entropy as a distance measure of cor- 
relations we investigate the important issue of evaluating quantum correlations such 
as entanglement, dissonance and classical correlations for 2”-dimensional Bell-diagonal 
states. We provide an analytical technique, which describes how we find the closest clas¬ 
sical states(CCS) and the closest separable states(CSS) for these states. Then analytical 
results are obtained for quantum discord of 2 n -dimensional Bell-diagonal states. As illus- 
tration, some special cases are examined. Finally, we investigate the additivity relation 
between the different correlations for the separable generalized bloch sphere states. 

Keywords: Quantum Discord, Distance Measure of Correlations, Dirae 7 matrices, 
Bipartite Quantum System. 

PACs Index: 03.67.-a, 03.65.Ta, 03.65.Ud 
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1 Introduction 

Quantum entanglement plays an important role in the quantum co mmuni cat,ion protocols 
like teleportation E EI, superdense coding [3], remote state preparation pt], cryptography 
[5] and many more. However, entanglement is not the only correlation that is useful for 
quantum information processing. Recently, it is found that many tasks, e.g. quantum non- 
locality without entanglement [SEI E], can be carried out with quantum correlations other 
than entanglement. It has been shown both theoretically and experimentally P ES that 
sorne separable states may speed up certain tasks over their classical counterparts. Recent 
measures of nonclassical correlations are motivated by different notions of classicality and op- 
erational means to quantify nonclassicality. One kind of nonlocal correlation called quantum 
discord, as introduced by Oliver and Zurek mm has received much attention recently 
[ISEIESESEIIESEIIESEIIESESIMIES. Most of these works are limited to studies 
of bipartite correlations only as the concept of discord, which rclics on the definition of mu- 
tual information, is not defined for multipartite systems. It is wcll known that the different 
measures of quantum correlation are not identical and conceptually different. For example, 
the discord does not coincide with entanglement and a direct comparison of two notions is 
rather meaningless. Therefore, a unified classification of correlations is in dernand. Modi et. 
al.|19j. introduced a unified classification of correlations for quantum states which is appli- 
cable for multipartite systems. Iu this unified view, the measure of correlation is based on 
the idea that the distance from a given state to the closest state without the desired prop- 
erty (e.g. entanglement or discord) is a measure of that property. Finding the CCS is stili 
a very difhcult problem and has the same challenged as faced in computing original discord. 
[261 27] . The examples of entangled states p with analytical expression for the CSSs, discussed 
in Refs.[28, [293, ED, EU [32, 33| Ei, EU SBj. The inverse problem to the long standing problem 
133 of finding the formula for the CSS was solved in [38JJ for the case of two qubits and a closed 
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formula for ali entangled States was solved in [39] for ali dimensions and for any number of 
parties. In this paper, we give an efficient procedure so that analytic evaluation of quantum 
discord of 2 n -dimensional Bcll-diagonal states can be performed. Then we find an exact ex- 
plicit formula for quantum discord of these states. We also show that total correlation for the 
separable generalized bloch spliere states is subadditive. 

This paper is organized as follows. In the next section we introduce the distance measure of 
correlation and show that the generic state p and its CSS or CCS have the sanie structure. 
In section 3 the dehnition of 2 n -dimensional Bcll-diagonal states is given and in section 4 we 
calculate the CCS of 2 n -dimensional Bcll-diagonal states and then we find an exact analytical 
formula for the quantum discord of these states. In the rest of this section to illustrate how 
the formula can be applied, we give two examples. In section 5 we investigate additivity rela- 
tions between different correlations for separable generalized bloch sphere states. Concluding 
remarks and two Appendices close this paper. 

2 Distance measure of correlations. 

Here, we will follow the approach of [T9j to characterize and quantify all kinds of correlations 
in a quantum state. The dehnitions of relevant quantities are: 


Entanglement E — min S'(/9 ||<j), 

a^T> 

(2.1) 

Discord D = min S(p||x), 
xec 

(2.2) 

Dissonance Q — min 1 S , (<j||x), 

(2.3) 

Classical correlations C = niinS'(x||7r), 

irGV 

(2.4) 


where V is the set of all product states (i.e., states of the form 7r = 7Ti ® 7t 2 <g)... <g) Kn and 7T n is 
the reduced state of the nth subsystem). C contains mixtures of locally distinguishable states 
X = Yjk n Pki...k N \ki...k N ){k 1 ...k N \ = J2kPk\k){k\ where p% is a joint probability distribution 
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and local states \k n ) span an orthonormal basis, V is the set of all separable states (i.e., states 
of the form a = EfcPfc 7r i ® ® ••• ® and 5'(a;||2/) = Tr(x logx — xlog y) is the relative 

entropy of x with respect to y. 

A bipartite state is called classical if it contains mixtures of locally distinguishable states, and 
is called separable if it can be represented as a convex combination of product states. Finding 
out the CSS is a non trivial task[2?j- While the set of separable states is apparently convex, 
this is not the case for the set of classical states then determining the CCS is even more 
complicated. Here we present an analytical procedure that allows us to obtain the CSS and 
CCS for 2 n -dimensional Bcll-diagonal states. The key idea is to hnd the minimum distance 
frorn a given state p to the set of all states without the desire property. The following theorem 
plays a Central role in minimizing the mentioned distance. 

Theorem: Given a generic state p G "H\X and A" G X, min S^pHX) is achieved when p and 
X have cornmon eigenbasis. Here, X is a special subset of the Hilbert space 'H. 

To show this, suppose 


N N 

p y ^ Aj|Aj)(Aj|, x y ^ Pj i py) (p,j i, 

* j 


then we have 


min5'(p||A) = min[Trplogp - Tr(plogX)] = A^ log A, - maxA^| (\j\pj ) [ 2 log pj (2.5) 


Suppose \(Xi\pj)\ 2 = qtj , where 


i,3 


= !, I> = 1- 

* 3 

Now the problem of finding the closest state A to p is reduced to the problem 

{ maximize Ey Aiqij log/ij = A T Qp. 
subject to E; Qij — Ej Tj = 1 
where p T = (logpi, log p 2 , log p N z), X T = (Ai, A 2 ,..., X n2 ) 


( 2 . 6 ) 


(2.7) 
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Eq. (12.61) shows that the matrix Q with the (Q)ij = Qij is doubly stochastic matrix. The 
set of doubly stochastic matrices, Q n , is the convex hull of the permutation matrices (Birkhoff 
(1946), von Neumann (1953)). In other words, the doubly stochastic matrix, Q n , is the convex 
combination of the permutation matrices, P n , that is 



( 2 . 8 ) 


so, Eq. (12. 7 j) takes form 



(2.9) 


maximize 


hence our problem reduces to a Linear Programming optimization over the convex set of feasi- 
ble region. Here the feasible region is a simplex and its apex, yicld when one of the r t equals 1 
and the others equal zero, are the desired Solutions of this optimization problem. This means 
that p and X have common eigenbasis. 

3 Definition of 2 n -dimensional Bell-diagonal states 

In order to put our discussion in a precise setting, let us first introduce 2 n -dimensional Bell- 
diagonal states acting on a bipartite systern H A ®H B with dim(H A ) = N = 2 n and dim(H B ) = 
N = 2 n . To do this, let S =< <71 = 71 < 8 > 71 , ...,<? 2 n = 72 n < 8 > 72 n} > be generated by 2n 
independent and commuting elcment such that -I is not an elcment of S and g'f = / for ali 
g l £S. 

'jj for j =1, 2, , 2n+l, known as Dirae matrices. (For a bricf review about Dirae matrices 

and an explicit construction of 7 j, see Appendix II). 

Hence, we can represent the density operators acting 011 a bipartite systern H A x H B as: 





(3.10) 
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where n,i 2 , e {0,1} 

Consicler the projectiori operators with 


2 n 


7 n 


U^ 2 v>* 2 - 


n 2 2n 




3 = 1 


Jlj2,-J2n 


,12, • • ■ ,*2n T?! J2, — J2n °iljl ''' ®i2n32n ,*2, — ,*2n 5 


E 


' K ii,i 2 ,...,i2n 


ll,l2,--.,*2n 


then we get 


p y ^ p*i,* 2 ,-") 


,*2n^”u,*2,---,i2n 


Jl,l2,...,l2n 


where 


(3.11) 


(3.12) 


■P*l,*2,--.,*2n 


ap E (- 1 )*” 


^ljl+^ 2 i 2 + ---+« 2 ni 2 n-/ 


v 3lj2> — j2n 


Jlj2,-J2 n 


T _ \ ^ /_l\iljl+i 2 i 2 +...+« 2 ni 2 n r) 

V jlj2, — d2n / j V - 1 -/ ni,i2r",i2n 

U)i2v^2n 

From the theorem above it follows that the CSS states can be represented as: 


(3.13) 


(7 = 


IV 2 


^ ] til,i2,—,i2n9l 92 ••■92n y ! Pil,i 2 ,-~, 


*2n^"*l,*2,.--!*2n 


*1,*2,...,*277 


*l)*2)-..)*2n 


(3.14) 


and CCS states are 


X P = Jp ^Ui2,:.,i2n9l 92 -92n = 


Il,l2,---,l2n 


.* 2 ,--.i* 2 n^"*l,* 2 ....,i 2 ri 

2l,*2,--.,*2n 


Xcr jy2 ^ 1 til,i2,---,i2n9l 92 ■■■92n 

*li*2.....*2n 


,l2,.--,*2n 7 hl,J2,...,*2n (3.15) 

11,12, —,*2n 


4 Calculation of classical states and quantum discord 

First of ali, we note that the resuit of above theorem can straightforwardly be used to obtain the 
CCS states of 2"-dimensional Bell-diagonal states. To do this recall that the Pauli operators 
on a single qubit are {/, cr x , <j y , a z j. The representation of the Pauli group we will deal with 
is the group forrned by elements of the forrn G n = { i k P\ <g) P 2 <8> ... <S> P n } where each P t is an 
element of {/, cr x , <j y , cr z j. 
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Suppose T is a subgroup of G n generated by elements {T i 1 i 2 ...i n jij 2 ...j n }- There is an ex- 
tremely useful way of presenting the generators of T [30]. To do this, we use r(rj li2 . Jn ) = 
\iii 2 --d n \jij 2 ---j n \ to denote the 2n-dimensional row vector representation of an elcment of the 
T. The left hand side of the row vector contains ls to indicate which generators contain o x s, 
and the right hand side contains ls to indicate which generators contain cr z s; the presence of a 
1 on both sides indicates a a y in the generator. More explicitly, it is constructed as follows. ff 
7 i contains an 1 on the jth qubit then the jth and n+jth column elements are 0; if it contains 
an cr x on the jth qubit then the jth column element is a 1 and the n+jth column element is a 
0; if it contains a o z on the jth qubit then the jth column element is 0 and the n+jth column 
element is 1; if it contains a a y on the jth qubit then both the jth and n+jth columns are 1. 
Let us define a 2 n x 2 n matrix A by 


0 In 


\ 


0 


\ 


/ 


(4.16) 


then the elements T 




i and T; ; < ■ ■ ■ are easily seen to commute if and only if 


r(T 


I\l2---ln]l32- 


= 0 ( m ° d 2 )' 


(4.17) 


Let g Ck G {r i 1 i 2 ...i n j 1 j 2 ...j n ® ^hi 2 -i n jij 2 ...j n } for k = 1, ...,n , such that the Eq. (|4.17l) is satisfied, 
then we can rewrite the state p such as: 


P = 


N 2 


E t n^ 1 n in n in+1 n^ 2n 

L n,---,i2 n yci •••y2n 




Since classical states contains mixtures of locally distinguishable states lience we can rewrite 


the expression of the CCS of p such as: 
1 


X = 


E 


jy2 / -> ^ii,i2,...,in,0, 0, 09ci"-9cn9n+l m --92n jy2 

nfold 


i 1 /-,in „0 


.0 




nfold 


(4.18) 


where = 'o_^o = i- 

nfold nfold 

Here we show that some of t- • n n n are zero and rest of them are the same of 


nfold 
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t i 1 12 in o o 0 - T° show this note that the eigenvalues of x which are 2 n —fold degener- 

nfold 

ate, are given by 


n _ _ \ ^ (_1 )hjl+-+i n jnf 

( Jtl,...,ln,in+l,---,l2n 7y2 / v V *~) ° A - A - 


Jl.-Jn 


31,32,...,3+0, 0, 0’ 

fold 


(4.19) 


so the problem of finding the CCS to p is reduced to the problem 

minimize S(p||y) 

SUbjeCt tO Eu^Wn+l,..,^ ^ 

The dual Lagrangian associated with this problem, is given by 


(4.20) 


^n+1 v>^2n 


L 'y 'j Pil,---,i2n^SPil,—,i2n ^ ^ ^ ^ Pii + l v,*2n) ^g + 1,... ,in,i-n 

Zi,...,Z 2 n 'i \5 - - - ^n+1 )***j^2n 

T p[ ^ 1 Qh,...,in,in+l, — ,i2n ~ 1] ' (^-21) 

ilr-^n^n+lr-i^n 

By calcnlating the gradient of the dual Lagrangian with respect to t jt ^ q q q and 


making it zero we get 
dL 


nfold 


dt ; 


- E 


x. 


n+l)-"^2n 


5^715^71-1-1 5’ ■ • 5^271 


31 J2v,jn, 0, 0,—,0 

nfold 


Q.i \ 5 ... 5 ^71 5^71+1 5 ••. 5 ^: 


+ ^](-l) iui+ - +w " = 0. (4.22) 


*2n 


Since 


E 


Pu,. •5^715^71+1 5' ■ * 5^271 


* I /fc 

jy2 zE ^ i Sl J2,..jn,0, 0, ..., 0’ 


^71+1 5'”5^2tI 


(4.23) 


3l,j2,---i32n 


n,fold 


then using Eq. (I4.22p one can show that 

tji,j 2 ,—,jn, 0 , 0,-.,0 — ^31,32,—,jw, 0, 0-j■ ■ ■ 0 ^ 1 ’ "">dm 

nfold nfold 

then we can rewrite the CCS of p such as: 


(4.24) 


X = 


1 h 

JX[2 / ./ ^* 1 , 12,..,in, 0, 0, ..., 0 + 1 - 


C 71 


(4.25) 


l\ 5 .. -5^271 


ifold 
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In the rest of this section we calculate the quantum discord of 2 2n -dimensional Bcll-diagonal 
states. Note that the key difference between the original definition of discord [TTj, 12] and the 
dehnition in Eq. fl2.2p is in minimization. We minimize the quantity D, while for the original 
discord, D — L p is minimized [E] where L p = S(n x ) — S(n p ) . Since for the Bcll-diagonal 
states L p = 0, hence the two forms of discord are the sanie. Now, using Eq. (l4.24p we give an 
exact analytical formula quantum discord for 2"-dimensional Bcll-diagonal states such as: 


D= Yl Ph,i 2 ,...,i 2 n l °gP 


u,*2,...,i 2 n - m ax — 


E E(-d 


iljl+—+injn 


*l,*2,.-..*2 n 


*l,*2,...,*n 31, -,3n 


E I” 1 ) 


iljl+---+injn 


tjl ,j2,-Jn,0, 0, .... 0 


1 , 32 , — ,jn , 0, 0, ..., 0 

nfold 

(4.26) 


3l,-,3n 


nfold 


men- 


where the maximum is taken over ali parameters ■ q 0 oi- Here we should 

nfold 

tion that the set {g l c \---g l c ™} in Eq. (14.251) can be chosen in many different ways. Since in the 
optimum strategy D is at its minimum, then we chose the set {< 7 **... < 7 ®"} or equivalently the 
parameters {tj 1J2 Jn () 0 o} suc h Hiat the second part in the Eq. (14.261) is maximized. To 

nfold 

give an intuitive understanding of this subject, let us illustrate it by a fundamenta! examples. 


4.1 Example 1: Generalized bloch sphere states 


Using Eq. fl3.10p the generalized bloch sphere states |4TJ are given by 

^ 2 n 

P= + 1 + tl,l,...,ig 2 n+l] 


(4.27) 


k =1 


where g 2n + 1 = 5 , i5 , 2---fl f 2n- Since classical states contains mixtures of locally distinguishable 
states hence using Eq. (14.251) one can conclude that the CCS of generalized bloch sphere states 
lic on the Cartesian axes. That is only one of the {t^ i2 iri 0 0 o) nonzer °- We assume 

nfold 

that, without loss of generality, U,o,...,o 0, then the CCS of generalized bloch sphere states 
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are given by 

X = jp[I + ti,o,...,o9i\, (4.28) 

In the optimum strategy -S*(p||x) is at its minimum, that is we have £i,o,...,o = t max - Hence 
using, (14.261) . we obtain 


N 2 — 1 

D = lo gPii,*a,...,i 2 n - --!og(l - tmax) ~ 1 + * nia * log(l + t max ) + 21og(JV). 

h,i2,-,i2n 

(4.29) 

This is in agreement with the resuit obtained in [4Tj for N=M. 

4.2 Example 2: 2 2 -dimensional Bell-diagonal states 

As the second example, to keep our discussion simple, let us focus attention on n =2 case. In 
this case, using (14.25P we have 

X = -^ E ‘i.,», 0 . 0 sZsS (4-30) 

where, using (14.171) we get 

(ScuSca) ^ {(riooo ® TiooO) r 0 ioi ® Toioi), (r 00 oi ® Toooii Tono ® Tono), (riooo ® riooo, Tono ® Tono), 
(r 0 ioo ® r 0 ioo, r 00 ii ® r 0 on), (rim <g) r 1111; r 1100 ® r 110 o), (rmi <s> r mi , r 10 io ® r 10 io), 

(riooo ® riooo, Toon ® Toon), (rooio ® Tooio, Toioi ® Toioi), 41111 <£) Tim, Tiooi < 8 > Tiooi), 

(Tooio ® Tooio, Tiooi < 8 > Tion), (roioo ® Toioo, Tioio Tioio), (roooi ® ^ 001 , Tnoo < 8 > Tnoo) 

(rooio ® r 0 oio,rnoo rn 0 o), (roioo ® roioo, riooi ® riooi), (roooi ® r 00 oi,ri 0 io <s> rioio)}, 

(4.31) 

and 

r 1000 = cr x < 8 > /, r 0 ioo — ® a xi Tooio — i& y ®cr y , r 0 ooi = i& y ® cr z ,r mi = —cr z <g) I, 

Tnoo = ~ (J z ® Ox, Tioio = — fz ® Tiooi = 2cr 2 <g) cr 2 , Tono) = —<g) a z , Toioi = I ® cr^, 

Toon = — 2 <g) da,, Toni = a x ® I, Tion = cr x { g> ctr, Tnoi = cr x ® a y , Tmo = cr x ® a z . 

(4.32) 
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Let us choose g ci = Tiooo riooo and g C2 = r 0 ioi <8> r 0 ioi, hence we have 


x — 77(1 + HoooTiooo ® r 1000 + ^oiooToioi ® r 0 ioi + tiioornoi ® r n0 i 
16 


Tlien, using (14.26[) quantum discord is given by 

4 

D = ^ ^ Ph,i2,i3M Pii, 12 , 13 ,-»4 ~ ^ ^ Pgi^°SPgi- 

*1,*2,*3>*4 i= 1 

where 

Pqi = illt + Hooo + £0100 + Hioo] 

= + ^1000 — hnoo ~~ ^1100] 

Pq 3 = [1 — Hooo + £0100 ~~ ^1100] 

Pq 4 = 1^[1 — HoOO — toioo + HlOo] 


(4.33) 


(4.34) 


(4.35) 


5 Subadditivity of correlations of a quantum state 


It has been conjectured [19] that the correlations of a quantum state are subadditive in the 
sense T p > E + Q + C c (where T p is total mutual information which dehned as S(p\\ti p ) 
and C a is the classica! correlation *S'(x 0 -|7r cr )). I 11 general, from an analytical point of view, 
the derivation of closed expressions of relative entropy of entanglement involves optimization 
procedures that are very complicated to perform. Hence, here we consider the inverse problem 
[38] and investigate additivity relations between different correlations for separable generalized 
bloch sphere states. These states are given by 


a = 


^ 2 n 

H' + £4. .1 ,...,0 9 k + il,l,...,l< 72 n+l]) 


N 2 


k =1 


(5.36) 


with the eigenvalues 


2 n 


^[i + E(-i)“4_ 


fc=i 


1 ,...,0 

k 




(5.37) 


where g 2 n+i = gig 2 ---g 2 n • The separable generalized bloch sphere states are actually bounded 
by Efc=1 . ol + l^i.i .il < 1 or > equivalently, p iu i 2 ,..., i2n < jp (Vii, i 2 ,..., i 2n ). The 
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family of all entangled states, p(x, a) , for which a is the CSS is given by [39] 


p(x, a) = a - xL a 1 (w iui2j ... ji2n ), 0 < x < x, 


max 


(5.38) 


where L a is linear operator. In the eigenbasis of a, a = diag(ai, ...,a n ) is a diagonal matrix, 
where ai,..., a n > 0 and for any /3 = [6” - =1 ], L a (j3) is defined by 



(5.39) 


u 'n,* 2 ,...,i 2 n; s are entanglement witnesses (EW) of 2 n -dimensional Bell-diagonal states. Here, 
Xmax is defined such that p(x max , cr) G H Uj+t i(convex set of positive hermitian matrices of trace 
one) and p(x max , a) has at least one zero eigenvalue. We also note that 'uy M2: ... : i 2 „. is normalized; 
i.e. Tr(iUi lt i 2} ,,, t i 2n ) 2 = 1 and Tr^L^ 1 = 0. In general EW of the 2 2n -dimensional 

Bell-diagonal states is given by [42] 



*l+*2 + ...+*2: 


lln+1 ® 7 2 2 n+l] 


(5.40) 


with the eigenvalue 



i 1 +^2 + • • • +^2n + j 1+J2 + • • • 


(5.41) 


Now dehne the real symmetric matrix [H^J 


[S( (T )]k 1 ...k2„,ll---h; 


Pk\...k2n Pll—hn. 


(5.42) 


1 n 4...fe„ - hl Ph...hJ 


lience we get 



(5.43) 


Note that L a is an invertible operator, where E 0 . 1 (tCj l! j 2v .. i j 2ri 
is the entrywise product of two matrices of A and B. Let 


) = • S(a) where A» B 



31—32 n 


(5.44) 
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then, using (15.38115.42115. 43 j) . we obtain 


Ph,j 2 ,.,hn = Ph,h ../,J - 


(5.45) 


We now focus our attention on the subadditivity of correlations. By direct calculation one gets 


E + Q + C a -T p = log ^ 1J2 ’- J2n + log l>n - 12 . Ja " 


Pjl,j 2 ,--,j 2 n 


( ljl,j 2 ,-,j 2 n 


+ 4jl,j2,-,j2n Qjl,j2,-,j2n y , Pjl,j2,—,j2n ^°g Pjl , J2, - - - J2n ^ y (Pjl ,32, -, J2r» Ph,32,- J2n ) 1°S Pjl,32, — J2n J 

(5.46) 


then, using (15.451) . we have 


E + C + CV-r,-* E dtSEZVn-n .+ . 


Jlj2v-J2 n 


Assume, with no loss of generality, the u>o,o,...,o> then we have that is 


2n 


y. 4. 1 — (-l) n ii,i,...,i + 1 — 0. 

t =i 


Thus, we have the following optimization problem 


max ,r ^ \tj 2 ):::yPii 42 ,..j 2 „logPnP 2 ,..,^ 


jl,j 2 ,-,j 2 n 


(5.47) 


(5.48) 


(5.49) 


subject to < 


3 

y=i4,...,^l^,...,o ~~ + 1 = 0 

3 

The dual Lagrangian associated with this problem, is given by 


2 n 


(5.50) 


L = X X h,S^Ph, 02 ,..., 32 n Ph,j 2 ,-..,j 2 n + M0,..., 0 [^ 4 ,..., 1 ,...,0 “ ( “ 1 ) ^1,1,..., 1 + 1] 

3 =1 


Jl,j 2 ,---,j 2 n 


2 n 

/^n,...,i 2ri+ i[y (~l) tj 4. 1 ,...,0 + ( _ 1)* 2 " +1 4,1,...,1 - 1], 

n,--.,22,1+1 j=l 


E 


(5.51) 
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and the complementary slackness condition(see the Appendix I) is given by 


2 n 


/*..wEHH.i.„+(-ir- + 'ii,i.i - 1 ] = o 

i =1 


(5.52) 


The possible optimal Solutions for above problcm are the edge and vertices Solutions. First of 
ali, we consider vertex solution that is one of the vertices. In the case of odd n, we have 


fi,i,...,i - -1 

4 ,..., 1 ,...,o = 0 > Vj G { 1 , 2nj 


(5.53) 


then Pj 1 ,j 2 ,-,j 2 n ~ 0 if ji + jz + ... + j-in— odd and Pj 1 ,j 2 ,...,j 2n — 772 if ji + J 2 + + j- 2 n— even. 


Hence Eq. (l5.47p gives 

E + Q + C a — T p = - 


2 n 


— (2n , 1),E log 1- Y [1 + Y(-l) jk - (— iyi+h+...+h n] = 0 
2 3 "-V 2 (n+l) N 2 . 1 f- v ' V ' J 

v v ' 31,32, — ,J2n k=l 

(5.54) 

Similarly, for other vertices( say = —1 ,t 0 ^ 0 = 0; Vj e {1,.... 2n}(even n) and 

3 

, with no loss of generality, t i,o,...,o = —1, — t Q ]_ 0 = 0) Vj G {2,..., 2 n}) one can 

3 

show that E + Q + C a — T p = 0. Let us next turn our attention to the edge Solutions. In this 


case we have 


L Pil,---,i2n^SPil,...,i2n 4 ^ Pil,---,i2n^SPil,—,i2n 

+ A t O,...,o[ 4 u )...0 + 4 ) 10 ...0 + + 4 ) 0 ...1 — (— l) n 4 u...l + 1 ] 


(5.55) 


By calculating the gradient of the dual Lagrangian with respect to making it zero one can 
show that 


Ho...o — hno...o — ••• — t, 


00 ... 1 - 


(5.56) 


Hence, using (15.3715.4815.56p . one finds that 

I 2n 

E + Q + C a — T p = 


2 3 V2(2» + o ( H 


£[i+ ^ ^ ^^ _j_ ^_j^H~Hi+...+t2n | 


h—fan A;=l 
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2 n _< 2 n 

(^(-l) ifc + 2 n(-iy i+h+ - +i Qt w ... 0 \ log ^[1 + (-1)' 1 '' 11 -"' 2 " + (^(-l) ifc + 2 n(-iy i+il+ - +in )t w ... 0 } 

k= 1 k =1 

(5.57) 

The above function is convex and it is zero at t 10 ...o — 0, then one can immediately deduce 
that E + Q + C c — T p becomes negative for the acceptable value of fio...o- In general, the 
Eq. (l5.47p is the difference of two convex functions, which both of thcm are non-positive and 
the intersection points of these two functions are in the vertices of the feasible region. On the 
other liand, for the edge and vertices Solutions E + Q + C a — T p < 0. Then one can conclude 
that the correlations of generalized bloch sphere states are subadditive in the feasible region. 

6 Conclusion 

In the unified view of quantum and classica! correlations the quantifications are done by the 
relative entropy and optimization of relative entropy is known to be a difficult problcm. In this 
work, we have presented a general algorithm via exact convex optimization to the problem of 
finding CSS and CCS for a given entanglement state p. Using the obtained CCS for the 2 n - 
dimensional Bell-diagonal states, we have derived analytical formula for the quantum discord 
of these states. As illustrating examples, we have analyzed the case of the separable generalized 
bloch sphere states and 2 2 -dimensional Bell-diagonal states and described how to apply the 
formula for this cases. We have also shown that the separable generalized bloch sphere states 
is subadditive. Wliile our analysis is for a special states of bipartite quantum system, it serves 
to provide a unified explanation for a variety of states . In fact, this approach is completely 
general and could be applicd for multipartite states in ali dimensions. The main conclusion 
is that the presented algorithm provide indispensable prerequisites for further investigation 
and can bring a robustness in constructing CSS and CCS for a given multipartite states in all 
dimensions. Application of this algorithm to other quantum system and finding related CSS 
and CCS is stili an open problem which is under investigation. 
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Appendix I: 

Convex optimization review: An optimization problem pfOj, has the Standard form 

{ maximize fn(x). 

(1-1) 

subject to fi(x) < bi,i — 1 ,m hi(x) — 0, i — 1, ...,p 

Where the vector x = (xi,...,x n ) is the optimization variable of the problem, the f u n et, ion 
/o : R n —» R is the objective fnnction, the fnnctions fi : R n —» R. i = 1 are the 

(ineqnality) constraint fnnctions, and the constants b\,b m are the lirnits, or bonnds, for the 
constraints. A convex optimization problem, is an optimization problem where the objective 
and the constraint fnnctions are convex fnnctions which means they satisfy ineqnality fi(ax + 
f3y ) < afi(x ) + Pfi(y ), for ali (x,y,a,/3) G R with a + /3 = l,a > 0, /3 > 0 and the equality 
constraint fnnctions lifx) = 0 mnst be affine (A set C G R n is affine if the line through 
any two distinet points in C lies in C). One can solve tliis convex optimization problem using 
Lagrangian duality. The basic idea in the Lagrangian duality is to take the constraints in convex 
optimization problem into account by augmenting the objective fnnction with a weighted sum 
of the constraint functions. The Lagrangian L : R n x R m x R p —> R associated with the 
problem is defined as 

m p 

L{x, A, v) = f 0 (x) + ^ Xifi(x) + ^ Vihi(x). (1-2) 

i=1 Z—1 

The Lagrange dual function g : R n x R m x R p —» R is defined as the minimum value of the 
Lagrangian over x: for A G R m , v G R p , 

g(X, v) — inf x&D L{x , A, u). (1-3) 

The dual function yiclds lower bonnds on the optimal value p* of the convex optimization 
problem, i.e for any A > 0 and any v we have 


g(X, v) <p*. 


(1-4) 
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The optimal value of the Lagrange dual problem, which we denote d*, is, by definitiori, the 
best lower bound on d * that can be obtained from the Lagrange dual function. In particular, 
we have the simple but important inequality 

d* < p*. (1-5) 


This property is called weak duality. If the equality d * = p * holds, i.e., the optimal duality 
gap is zero, then we say that strong duality holds. If strong duality holds and a dual optimal 
solution (A*, v*) exists, then any prinral optimal point is also a mininrizer of L(x, A*, u*). This 
fact sonretimes allows us to compute a primal optimal solution from a dual optimal solution. 
For the best lower bound that can be obtained from the Lagrange dual function one can solve 
the following optimization problem 


{ maximize g(\,v). 
subject to A > 0 


( 1 - 6 ) 


This problem is called the Lagrange dual problem associated with the main problem. Condi- 
tions for the optimality of a convex problem is called Karush-Kuhn-Tucker (KKT) conditions. 
If fi are convex and hi are affine, and x, A, fi are any points that satisfy the KKT conditions 


hfx) = 0, i = 1 , ...,p, 

fi(x) <0, i = 1,..., m, 

A* > 0 A ifi(x) = 0, i = 1,..., m, 

m p 

V fo(x) + Ai V fi(x) + "i V hi(x) = 0. (1-7) 

i i 

then x and (A, fi) are primal and dual optimal, with zero duality gap. In other words, for 
any convex optimization problem with differentiable objective and constraint functions, any 
points that satisfy the KKT conditions are primal and dual optimal, and have zero duality 
gap. Hence, f 0 (x) = g(X, fi). 
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The condition A ifi(x) = 0, i — is known as complementary slackness; it holds for any 

primal optimal x and any dnal optimal (A, v) (when strong duality holds) 


Appendix II: 

Throughout the paper, we have used the formalism of Dirae 7 matrices. Therefore, in this 
appendix we dehne the algebra of Dirae 7 matrices and exhibit matrices wliich realize the 
algebra in the Euclidean representation and explain onr notations and conventions. 

To do this, let 7 /i = 1, d, be a set of d matrices satisfying the anticommuting relations: 

7^7 v + Ivi» = 2 5^1, ( 1 - 1 ) 

in which 1 is the identity matrix. These matrices are the generatores of a Clifford algebra 
similar to the algebra of operators acting 011 Grassmann algebras. It follows frorn relations 
(II-II) that the 7 matrices generate an algebra which, as a vector space, has a dimension 2 d . 
In the following, we will give an inductive construction (d —» d + 2) of hermitian matrices 
satisfying (II- II) . In the algebra one element plays a special role, the product of ali 7 matrices. 
The matrix 7 ,,: 

7s = *^7i72---7d> (1-2) 

anticommutes, because d is even, witli all other 7 matrices and 7^ = I. 

In calculations involving 7 matrices, it is not always necessary to distinguish frorn other 
7 matrices. Identifying thus 7* with 7 ^+ 1 , we have: 

lilj + Ijli = = 1,..., d,d+ 1. (1-3) 

The Greek letters /rz/... are usually used to indicate that the value d+l for the index has been 
excluded. 

An explicit construction of 7> 
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It is sometimes useful to have an explicit realization of the algebra of 7 matrices. For d — 2, 
the Standard Pauli matrices realize the algebra: 


(d= 2 ) 

7i =01 = 


1 0 

if 

(d= 2) ^ 1 


V 1 


,72 = (72 = 

* 0 ) 


(d=2) (d=2) 

7s - 73 = ^3 = 


^ 1 0 ^ 


V 


0 


(1-4) 




The three matrices are hermitian, i.e., 7 \ = 7 J. The matrices 71 and 73 are symmetric and 72 
is antisymmetric, i.e., 71 = 7 *, 73 = 73 and 72 = — 72 - To construet the matrices for higher 
even dimensions, we then proceed by induction, setting: 


(d+2) _ (d) 

7i = <71 ® 7i = 


( 0 7i (d) ^ 


V 


# 0 


, f — 1 , 2 ,d + 1 , 


/ 


7rf+2 = 0-2 ® I {d} = 


( n ■ T \ 

0 -il d 


V 


Hd 0 


(1-5) 


/ 


where, Id is the unit matrix in 2 2 dimensions. As a consequence 7 i d+2) has the form: 

\ 

7 ^ +2) = ld+3 ] =°3®h = 


( 1 - 6 ) 


/ 


-fd 0 
0 —Id 

A straightforward calculation shows that if the matrices 7 ^ satisfy relations (11-311 . the ”f\ d+2 ^ 
matrices satisfy the same relations. By induction we see that the 7 matrices are all hermitian. 
from (11-51) . it is seen that, if 7 ^ is symmetric or antisymmetric, has the same property. 

The matrix ^+ 2 ^ i s antisymmetric and the matrix 7^3 is symmetric. It follows immediately 
that, in this representation, all 7 matrices with odd index are symmetric and all matrices with 
even index are antisymmetric, i.e., 


7- = (-l) i+1 7n 


(1-7) 
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